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Abstract -We investigate, using Monte Carlo simulations, the phase diagram of a system of hard 
rectangles of size m x mk on a square lattice when the aspect ratio k is a non-integer. The existence 
of a disordered isotropic phase, a nematic with only orientational order, a columnar phase with 
orientational and partial translational order, and a high density phase with no orientational order 
is shown. The high density phase is a solid-like sublattice phase only if the length and width of the 
rectangles are not mutually prime, else, it is an isotropic phase. The minimum value of k beyond 
which the nematic and columnar phases exist are determined for m = 2 and 3. The nature of 
the transitions between different phases is determined, and the critical exponents are numerically 
obtained for the continuous transitions. 
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Understanding the nature of the different phases and 
the transitions between them in a system of hard rods has 
significance for more complex physical systems such as liq¬ 
uid crystals [l], tobacco mosaic virus [2], fd virus m, 
silica colloids EEj, boehmite particles 019], DNA origami 
nanoneedles 1101, and adsorption of gas particles on metal 
surfaces mm- In three dimensional continuum, the sys¬ 
tem of hard rods undergoes an entropy driven phase tran¬ 
sition from a low density isotropic phase to a high den¬ 
sity nematic phase that has orientational order 1X741X8] . 
Further increase in density may result in a smectic phase 
with orientational and partial translational order and a 
solid phase [19] . In two dimensions, the system under¬ 
goes a Kosterlitz Thouless type phase transition from an 
isotropic phase to a power law correlated phase [20]. Hard 
cuboids on a cubic lattice and hard rectangles on a square 
lattice are the corresponding lattice analogues. While the 
complete phase diagram for cuboids is not known, the sys¬ 
tem of hard rectangles has a rich phase diagram m- 
Consider a system of monodispersed hard rectangles of 
size m x mk (fc > 1) on a square lattice, interacting only 
through excluded volume interaction. When to = 1 (hard 
rods), the system undergoes two transitions with increas- 

( a )joyjit@imsc.res.in 

T) rrajesh.imsc.res.in 


ing density for aspect ratio k > 7: first from a low den¬ 
sity disordered phase to an intermediate density nematic 
phase [22H261 . and second, from the nematic phase to a 
high density disordered phase [27][28]. The high density 
disordered phase has been argued to be a reentrant low 
density disordered phase [29]. When k = 2 (dimers), it 
may be rigorously shown that there are no transitions 
at non-zero densities of vacancies 137711551 . though at full 
packing, the correlations decay algebraically on bipar¬ 
tite lattices 134) . and exponentially on non-bipartite lat¬ 
tices [35]. For k 1, the existence of the nematic phase 
at intermediate densities may be proved rigorously [231 . 
The only exact solutions that exist are for the model of 
hard rods on a Bethe-like lattice [29ll36l . 

For m > 2 and integer k, four different phases have been 
observed: isotropic, nematic, columnar, and solid-like sub¬ 
lattice phases [2X1571155] . The nematic phase exists only for 
k > 7 for to = 2,3. The columnar phase exists only when 
k > 4 for to = 2 and for k > 2 for to = 3. For large enough 
k , with increasing density, the system transits successively 
from isotropic to nematic to columnar to sublattice phase. 
The nature of the phase transitions has also been studied. 
The isotropic-nematic transition belongs to the Ising uni¬ 
versality class for all to. When to > 3, all other transitions 
are first order in nature. When to = 2, all transitions ex- 
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cept the isotropic-columnar transition are continuous and 
belong to the Ising or Ashkin-Teller universality class. The 
isotropic-columnar transition is continuous for k = 5 and 
first order for k = 6, implying the existence of a tricritical 
point at an intermediate value of k [2T :. 

All the above results are for integer values of the as¬ 
pect ratio k. Within viral expansion and Bethe approx¬ 
imation, the results do not depend on whether k is an 
integer or not m ■ However, these results are presum¬ 
ably valid only for large k. For smaller k , when the ap¬ 
plicability of the model for adsorption of gas particles on 
metal surfaces is more relevant, it is not clear what the 
effect of k being an non-integer is. What are the different 
phases and the phase diagram when k is rational but not 
an integer? What are the minimum values of k beyond 
which the nematic and columnar phases exist? Answer¬ 
ing these questions will allow us to obtain the complete 
phase diagram for the system of hard rectangles. In this 
paper, we obtain the phase diagram for m = 2, when k is 
a half-integer, using large scale Monte Carlo simulations. 
The isotropic-columnar transition is shown to be discon¬ 
tinuous while the isotropic-nematic and nematic-columnar 
transitions are shown to be continuous. The critical expo¬ 
nents are numerically estimated for the continuous transi¬ 
tions. We find that the columnar phase exists only when 
k > 11/2 for m = 2 and when k > 13/3 for m = 3. The 
nematic phase exists only when k > 15/2 for to = 2 and 
when k > 22/3 for m = 3. 

The Model and Monte Carlo algorithm. — Con¬ 
sider monodispersed hard rectangles of size to x mk on a 
square lattice of size L x L with periodic boundary con¬ 
ditions. Each rectangle is oriented either horizontally or 
vertically. A horizontal (vertical) rectangle occupies mk 
lattice sites along x (y)-direction and m lattice site along 
y (x)-direction. Each site may have at most one rectangle 
passing through it. We associate an activity z = e 11 to 
each rectangle, where p is the chemical potential. In this 
paper, we restrict the aspect ratio k to non-integers. 

We simulate the system in the constant p grand canoni¬ 
cal ensemble using an efficient algorithm involving cluster 
moves that has been shown to be very useful in equili¬ 
brating hard core systems of extended particles at high 
densities EUEziEg. Here we briefly review the algo¬ 
rithm. Starting from a valid configuration, a row or a 
column (say a row) is chosen at random. All horizon¬ 
tal rectangles whose bottom-left corners (heads) are on 
that row are evaporated, keeping the rest of the config¬ 
uration unchanged. The row now consists of intervals of 
empty sites, separated by sites that are either occupied 
by rectangles or can not be occupied due to the hard 
core constraint. The empty intervals of the row are re¬ 
occupied by a new configuration of horizontal rectangles 
with the correct equilibrium grand canonical probabilities. 
The calculation of these probabilities reduces to a solvable 
one-dimensional problem. If a column is chosen, similar 
evaporation-deposition moves are performed for vertical 


rectangles. Equilibration is faster on including a flip move 
in which a square plaquette consisting of £ aligned hori¬ 
zontal or vertical rectangles is rotated by 7 t/ 2, where l is 
the ratio of the least common multiple of m and mk to 
to. A detailed description of the implementation of the 
algorithm for the system of hard rectangles is described 
in Ref. EH- Unlike for integer k , the flip move is less 
effective for non-integer k because the rotatable plaque- 
ttes are larger in size and hence, have lower probability 
to occur during the simulations. This makes it difficult 
to equilibrate the systems at high densities. Other imple¬ 
mentations of the algorithm include lattice models of hard 
rods E3HH1, hard discs [35], and mixtures of dimers and 
hard squares [40] . 

Different Phases. — As for integer fc, we observe four 
different phases in the simulations: an isotropic (I) phase, 
a nematic (N) phase, a columnar (C) phase and a high 
density (HD) phase. The I phase is disordered. In the 
N phase, rectangles orient preferably along the horizontal 
or vertical direction, but they do not have any positional 
order. Each row or column on an average contains equal 
number of heads of rectangles. The columnar (C) phase 
has orientational order and translational order only in the 
direction perpendicular to the nematic orientation. When 
to = 2, in the columnar phase, if the majority of the rect¬ 
angles are horizontal (vertical), their heads lie mostly on 
even rows (columns) or odd rows (columns). Hence, there 
are 4 (in general 2 to) symmetric C phases. The I, N and 
C phases are observed when m > 2, for both integer and 
non-integer k. 

The HD phase has no orientational order. But it may 
or may not possess translational order depending on the 
length and width of the rectangles. Let, the greatest com¬ 
mon divisor of the length and width be denoted by p. We 
divide the square lattice into p 2 sublattices by assigning 
to a site (*, j) a label (* mod p) + p x (j mod p). In the 
fully packed limit, it is straightforward to verify that the 
heads of the rectangles occupy one of the p 2 sublattices. 
We expect this phase to be stable to introduction of va¬ 
cancies at densities close to the full packing. If p > 1, the 
HD phase is a sublattice phase with complete translational 
order but no orientational order [21j . On the other hand, 
when p = 1 (length and width are mutually prime), the 
HD phase is disordered with no orientational or transla¬ 
tional order. Since existing evidence for m — 1 suggests 
that the high density disordered phase is qualitatively sim¬ 
ilar to the low density I phase [281129] , we expect the same 
to hold for to > 2 whenever the HD phase is disordered. 

When to > 2 and integer k, p = m > 1 and the HD is 
known to be a sublattice phase [2Tj . consistent with the 
above argument. To further confirm that the HD phase 
is a sublattice phase when p > 1, but k is a non-integer, 
we simulate the system of rectangles of size 4x6, for 
which p = 2. We divide the lattice into p 2 = 4 sub¬ 
lattices. The sublattice order parameter is defined as 
qi = no — n\ — ri 2 + ri 3 , where tii is the fraction of sites oc- 
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Fig. 1: Time evolution of the sublattice order parameter gi, 
starting from three different types of initial configurations: ne¬ 
matic, disordered, and sublattice phases. The straight line is 
|qi| = 0.580. The data are p = 10.8 and for L=960. The 
equilibrium density is « 0.962. 

cupied by rectangles whose heads are on the ith sublattice. 
It is straightforward to check that (|(?i|) 7 ^ 0 only for the 
sublattice phase. To show the existence of sublattice phase 
at high density, a large value of p is chosen ( p = 10 . 8 ), 
and the temporal evolution of \q\ | is tracked, starting from 
three different initial configurations: nematic, disordered 
and sublattice phases (see Fig [T]). At large times, the sys¬ 
tem reaches a stationary state that is independent of the 
initial configuration, ensuring equilibrium. For this choice 
of /x, the fraction of occupied sites p fluctuates around 
0.962. In equilibrium, (|qi|) « 0.580, clearly showing the 
existence of a sublattice phase. 

Phase Diagram for m = 2. — The phase diagram 
for m = 2 and non-integer k is shown in Fig. [2j where the 
data points are obtained from Monte Carlo simulations 
and the lines, based on analysis of the phase diagram for 
large k m are guides to the eye. The low density phase 
is an I phase for all k. Since the length and width of the 
rectangles are mutually prime, p = 1, and the HD phase 
is a reentrant I phase. No phase transitions are observed 
when k < 9/2. The C phase exists only for k > 11/2, 
while the N phase exists only for k > 15/2. 

We could not numerically obtain any data point on the 
C-HD phase boundary as it is not possible to equilibrate 
the systems within available computer time at high densi¬ 
ties for k > 11/2. However, the critical density for C-HD 
phase transition was argued to behave asymptotically as 
1 — a/(mk 2 ), for k 1, where a is a constant. Likewise, 
for large fc, the critical density for the I-N phase transition 
scales as Afc -1 , where A is a constant, independent of m, 
and that for the N-C phase transition tends to a non-zero 
constant mm- The solid lines in Fig. [2] follow these 
asymptotic behavior for large k. 

The I-C transition is found to be first order for both 
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Fig. 2: Phase diagram for rectangles of size 2x2 k, where k is re¬ 
stricted to non-integer values. I, N, C and HD denote isotropic, 
nematic, columnar and high density phases respectively. The 
data points are from simulation, while the continuous lines and 
shaded portions are guides to the eye. The shaded portion de¬ 
notes regions of phase coexistence. 

k = 11/2 and 13/2. The shaded region in Fig. [2] denotes 
the region of phase coexistence at a first order phase tran¬ 
sition. We find that the I-N and N-C transitions are both 
continuous. These transitions are analyzed in detail be¬ 
low. 

Critical behavior for m = 2. — We now study the 
nature of the different phase transitions for the system of 
2 x 2k rectangles, where k is half-integer. To study the 
I-N transition, we define the order parameter 

q 2 = n h ~ n v , ( 1 ) 

where rih and n v are the fraction of sites occupied by the 
horizontal and vertical rectangles respectively. In the I 
phase (|< 72 1) = 0, while in the N phase, (|< 72 1) 7 ^ 0. 

The I-C and N-C phase transitions are best studied with 
the order parameter 

93 = | ^re Hro\ \^ce ^coli (2) 

where n re ( n ro ) is the fraction of sites occupied by rectan¬ 
gles whose heads are in the even (odd) rows, and n ce ( n co ) 
is the fraction of sites occupied by rectangles whose heads 
are in the even (odd) columns. In the I and N phases, 
n re ~ n ro , and n ce ~ n co , implying that (| 1) = 0. In the 

C phase, either n re 7 ^ n ro and n ce ~ n co , or n ce 7 ^ n co and 
n re ~ n ro implying that (I 93 I) 7 ^ 0 . 

The other relevant thermodynamic quantities are the 
second moment Xi and the Binder cumulant Ui , defined 
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as 


Ui 

Xi 


= 1 - 


(qf) 

3<<? 2 } 2 ' 

k 2 )^ 2 , 


(3a) 

(3b) 


where i = 2,3. Near the critical point, the singular be¬ 
havior is captured by finite-size scaling: 


<N> 


(4a) 

Ui 

^ fu{eL 1/u ), 

(4b) 

Xi 


(4c) 


where e = (p — p c )/p c , where p c is the critical chemical, 
(3, 7 , v are the critical exponents, and f q , /„, and f x are 
scaling functions. 

Isotropic-Nematic (I-N) transition. We study the I- 
N transition for 2 x 15 (k = 15/2) rectangles using the 
order parameter g 2 . Since the N phase may have orien¬ 
tational order only in the horizontal or vertical direction, 
we expect the I-N transition to be in the two-dimensional 
Ising universality class, as has been confirmed for integer 
k, when m = 1 [24] and m = 2,3 EB, and for systems of 
polydispersed rods mm- The data for U 2 for different 
system sizes intersect at p = pj_ N « 0.945 [see Fig.EJa)]. 
The corresponding critical density is pj_ N « 0.694, which 
is less than Pj_ N ~ 0.745 for k = 7 [28] , consistent with 
Pi_ N ~ Ak~ x [37] - The data for U 2 [see Fig. E[b)], (|g 2 |) 
[see Fig. He)], and \2 [see Fig. Eld)] for different system 
sizes collapse onto a single curve when scaled as in Eq. m 
with Ising exponents ft/v = 1/8, 7 / 1 / = 7/4, and v = 1. 
For larger values of k, integer or otherwise, we expect the 
I-N transition to be in the Ising universality class. 

Nematic-Columnar (N-C) transition. We study the 
N-C transition for rectangles of size 2 x 15 using the or¬ 
der parameter 53 . When the system makes a transition 
from the N phase with horizontal (vertical) orientation to 
the C phase, the symmetry between even and odd rows 
(columns) is broken. From symmetry considerations, we 
expect the N-C transition to be in the Ising universality 
class. The data for t /3 for different system sizes intersect 
at p = p c N _ c ~ 1.696 [see Fig. E[a)], corresponding to 
Pn~c ~ 0.759. The data for t /3 [see Fig. EJb)], (| < 73 1) [see 
Fig- Sic)], and y '3 [see Fig. Eld)] for different system sizes 
collapse onto a single curve when scaled as in Eq. (j4j with 
Ising exponents ft/v = 1/8, 7 / 1 / = 7/4, and v = 1. The 
N-C transition in the system of 2 x 14 rectangles has also 
been shown to be in the Ising universality class m , and 
we expect the same for k > 15/2. 

Isotropic-Columnar (I-C) transition. The I-C transi¬ 
tion occurs only for rectangles of size 2 x 11 (k = 11 / 2 ) 
and 2 x 13 (k = 13/2). The transition is best studied us¬ 
ing the order parameter q 3 . We find that the transition is 
first order for both values of k. This may be established by 
numerically calculating the probability density functions 
(pdf) P(p) of the density p and P{qz) of the order pa¬ 
rameter (73 for values of p that are close to the transition 
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Fig. 3: The critical behavior near the I-N transition for rect¬ 
angles of size 2 x 15 (fc = 15/2). (a) The data for Binder 

cumulant for different system sizes intersect at fs 0.945 

( Pi-n ~ 0 .694). The data for different L near the I-N tran¬ 
sition for (b) Binder cumulant, (c) order parameter, and (d) 
second moment of the order parameter collapse onto a sin¬ 
gle curve when scaled as in Eq. a with the Ising exponents 
fi/v = 1/8, 7 / 1 / = 7/4, and v = 1. 
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Fig. 4: The critical behavior near the N-C transition for rect¬ 
angles of size 2 x 15 (k = 15/2). (a) The data for Binder 

cumulant for different system sizes intersect at Pn-c ~ 1-696 
( Pn-c ~ 0.75). The data for different L near the N-C tran¬ 
sition for (b) Binder cumulant, (c) order parameter, and (d) 
second moment of the order parameter collapse onto a sin¬ 
gle curve when scaled as in Eq. a with the Ising exponents 
fi/v = 1/8, x/v = 7/4, and v = 1. 
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Fig. 5: Probability density function of (a) density p and (b) 
order parameter qs for three values of p near the I-C transition. 
The data are for rectangles of size 2x11 and L = 440. 

point pf_ c « 3.885 for k = 11/2 and p, ~ pj_ c = 2.390 
for k = 13/2. The pdf for k = 11/2 and 13/2 are shown 
in Fig. [5] and Fig. [ 6 ] respectively. In both the figures, the 
pdfs for have three clear peaks at the transition point: 
the two peaks at q% ^ 0 correspond to the symmetric C 
phases and the one at <73 = 0 to the I phase. The pdf for 
p have two peaks of nearly equal height at the transition 
point, though these peaks are not clearly separated for 
k = 13/2. The difference in the peak positions is equal to 
the jump in the density across the transition point and is 
shown by the shaded region in Fig. [2] We find that these 
peaks become sharper with increasing system size. These 
are clear signatures of a first order transition. 

Phase Diagram for m > 3. — We expect that the 
phase diagram for m > 3 and non-integer k to be qual¬ 
itatively similar to that for m = 2 with three entropy 
driven transitions for large k. The HD phase is a disor¬ 
dered or sublattice phase depending on whether the length 
and width of the rectangles are mutually prime or not. For 
m = 3, we determine the the minimum value of k beyond 
which the C and N phases exist. There are no transitions 
for k < 11/3. We find that the C phase exists only for 
k > 13/3, while the N phase exists only for k > 22/3. 
When m = 3, the C phase has a 6 -fold symmetry, and 
thus with analogy with Potts model, we expect a first or¬ 
der transition. For rectangles with k = 13/3, we confirm 
that the I-C phase transition is first order in nature. The 
pdfs for p and <73 behave similarly to that for the case of 
m = 2. All the other transitions are also expected to be 
first order except the I-N transition, as seen for integer 

k mj. 

Summary and Discussions. — In this paper, we ob¬ 
tained numerically the phase diagram of the system of 
hard rectangles of size m x mk with non-integer aspect 
ratio k. As for integer fc, the system may exist in four dif¬ 
ferent phases: isotropic, nematic, columnar or high density 
phase. For integer k , the high density phase is a solid-like 
sublattice phase. However, when k is a non-integer, the 
high density phase is a disordered phase when the length 
and width of the rectangles are mutually prime. The phase 
diagram for large m is expected to be qualitatively similar 
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Fig. 6: Probability density function of (a) density p and (b) 
order parameter q$ for three values of p near the I-C transition. 
The data are for rectangles of size 2 x 13 and L = 416. 

to that for m = 2. The isotropic-nematic transition will 
be in the Ising universality class for all m. However, all 
other transitions are expected to be first order. 

The N phase is found to exist only when k > 15/2 for 
m = 2 and when k > 22/3 for m = 3. For integer k , the N 
phase exists only when k > 7 for m = 1.2.3 I2T1I221 . These 
different lower bounds may be combined to give tighter 
bounds for the smallest value of k beyond which 

the N phase exists. We conclude that 20/3 < k^~^ < 7. 

The bounds for the minimum value of k beyond 

which the C phase exists, are not so clear. We find that 
the C phase exists when k > 11/2 for m = 2 and when 
k > 13/3 for m = 3. On the other hand, for integer k : 
the C phase exists for k > 4 for m = 2 and when k > 2 
for m = 3. Thus, unlike for the N phase, k^~^ depends 
both on m and whether k is a integer or not, and it is not 
possible to combine the bounds only in terms of m. 

The phase diagram when the aspect ratio of the rectan¬ 
gles is irrational remains an open question. In this case, 
it has been conjectured that there could be more tran¬ 
sitions at densities close to full packing, when the disor¬ 
dered phase will become unstable to a nematic or colum¬ 
nar phases [22]. This question, as well as finding tighter 
bounds on k for existence of different phases are best an¬ 
swered by obtaining the phase diagram of hard rectangles 
with restricted orientation in two dimensional continuum, 
when the aspect ratio may be continuously tuned. 

* * * 

The simulations were carried out on the supercomput¬ 
ing machine Annapurna at The Institute of Mathematical 
Sciences. 
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